
function Results = print_results_theta(results,dnames,fe)

KR = length(fe);
KD = length(dnames);

colNames = {'Estimate', 'StdError', 'pValue'};

rowNames = [cellfun(@(t)['PVEM x ',t],dnames,'UniformOutput',0),...
    cellfun(@(t)['PRI x ',t],dnames,'UniformOutput',0)];
rowNames = [{'prob_winning'},rowNames];
results([1:KR,KR+KD+(1:KR)],:) = []; % remove FEs
results = results([end,1:end-1],:); % rearrange

Results = array2table(results,'RowNames',rowNames,'VariableNames',colNames(1:size(results,2)));
Results.Properties.DimensionNames{1} = 'Coefficient';



